% return matrix (val) of size (siz) with i.i.d. uniform variables in [a,b]
function val = DrawUniform(siz,a,b);
val = rand(siz);
val = val*(b-a)+a;

